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Abstract. Using arguments parallel to those used in support of using H2CO 
as a sensitive probe of temperature and density in molecular clouds, we measured 
the J=7^6 and J=10^9 transitions of thioformaldehyde (H 2 CS) in several hot 
core sources. The goal here was to investigate more closely the conditions giv- 
ing rise to H2CS emission in cloud cores containing young stars by modelling 
several transitions. The H 2 CS molecule is a slightly asymmetric rotor, a heavier 
analogue to H2CO. As in H2CO, transitions occur closely spaced in frequency, 
though they are substantially separated in energy. Transitions of H 2 CS originat- 
ing from the K=0, 1, 2, 3, and 4 ladders in the 230 and 345 GHz windows can 
productively be used to constrain densities and temperatures. As a first step in 
developing the use of these transitions as thermometers and densitometers, we 
surveyed and modeled the emission from well known warm dense cores. 



1. Introduction 



The H2CS molecule is asymmetric and planar as is the case for H2CO; it has C2 V 
symmetry. The molecular asymmetry results in a substantial dipole moment 
of 1.65 Debye. An advantage of the spectrum of a slightly asymmetric rotor 
molecule is that lines arising between levels of different energies lie nearby in 
the spectrum, affording a gauge of excitation while observing with a single or 
only slightly different tuning on one telescope during one session. Calibration 
uncertainties are thereby minimized. In a single spectral region the lines of both 
formaldehyde and thioformaldehyde occur between two states with the same 
angular momentum quantum number, J. As there may be numerous values for 
the projection of the angular momentum vector on the main axis of the molecule, 
there are many rotational states with a single angular momentum quantum 
number J; the quantum number denoting these projections are usually denoted 
K. Except for the case that K=0, the spins of the hydrogen atoms further split 
each level denoted by a particular value of K into two levels. The spins of the 
hydrogen atoms also segregate two states of the molecule, ortho states in which 
the value of K is odd, and para states in which the value of K is even. As a result 
there are often an abundance of lines in any broad frequency window which arise 
from many states of excitation. Measuring these lines can result in a measure 
of the state of excitation of the gas in which they arise, in a similar manner to 
that for which H2CO lines are used (Mangum and Wootten 1993). 
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Figure 1. DR21(0H) spectra of the H 2 CS lines at 342.9 GHz (K=0 center) , 
343.8 GHz and 343.3 (K=2, upper and lower respectively of the bottom lines); 
and 343.5 GHz (K=l, upper) and 343.4 GHz (K=3, second from top; the left 
vertical flag marks the location of the second component of the blend) . The 
right vertical flag marks the adopted frequencies (from Splatalogue/SLAIM). 



2. Observations 

Measurements of a variety of transitions in the 230 and 345 GHz atmospheric 
windows were made on 2006/04/29-05/04 using the facility receivers of the 
Caltech Submillimeter Observatory (CSO). Pointing was checked several times 
each night and was found to be stable to within 5". Antenna efficiency was 
checked using position-switched measurements of Jupiter (O — 44.6, Tjupiter ~ 
170 K). At 240 and 345 GHz the CSO main beam FWHM measures ~30"and 
~21"respectively. Using the formalism described in Mangum (1993) we derive 
ilmb ~ 0.71 and ~0.61 at 240 and 345 GHz, respectively; results are given in Ta- 
ble 1. Lines observed in DR21(OH) are illustrated (Fig. 1); line ratios measured 
toward the source and derived parameters are given in Table 2. 

The lines in Table 1 cluster in frequency, allowing a single suitable band- 
width receiver to measure them all simultaneously or with little retuning (hence 
with a more accurate calibration). However, their energies are substantially dif- 
ferent. In one observation per band, then, one can observe several line ratios. 
As H2CS has a lower abundance than H2CO, and because there are more energy 
levels in any given interval, the molecular populations are spread thinly and the 
emission is expected to be optically thin. Line ratios should, therefore, provide 
a fairly robust estimate of excitation. For most sources in our survey, we have 
observed several lines at 1.3mm and several additional lines at .85mm. For a 
convenient excitation measure we employ line ratios, for the 1.3mm para-H2CS 
lines we define RqJ as the 7o7~6o6 to 726,5 ~ 625,4 line ratio, given in Table 2 
(where we have averaged the K=2 lines); a similar quantity for the .85mm para- 
H 2 CS lines, R^ 2 10 = T^(10 ,io - 9o,g)/ T^(10 2 , 8 - 9 2 , 7 ), where again the K=2 
lines are averag ed, and for the 0.85mm ortho-H 2 CS lines R/ 3 10 = T^(10i 9 - 

9i )8 )/ lfl(10 3 ,7 - 9 3 , 6 ). 
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Tabic 1. H 2 CS Results: DR21(0H) J2000 20:39:01.1 42:22:48 



Transition 


Frequency 


T* A 


V LSR 


AV 


R 




(GHz) 


(K) 


(km/s) 


(km/s) 


(K) 


10o,io — ► 9o,g 


342.946456 


0.32 ±0.01 


-4.5 


5.8 


0.30 


IO28 —> 927 


343.321887 


0.24 ±0.02 


-4.4 


6.1 


0.07 


l°3,8/7 9,3,7/6 


343.413647 


0.17 ±0.01 


-4.4 


5.2 




IO29 - ► 928 


343.812941 


0.33 ±0.02 


-5.1 


6.9 


0.07 


10ig 9i 8 


348.53424 


0.58 ±0.10 


-5.0 


6.5 


0.61 


7o7 -»■ 606 


240.26632 


0.49 ±0.02 


-3.2 


5.9 


0.51 


726 625 


240.38175 


0.12 ±0.03 


-4.0 


7.1 


0.13 


73,5/4 — *• 63,4/3 


240.39281 


0.15 ±0.03 


-4.0 


5.0 


0.11 


725 — ► 624 


240.39348 


0.13 ±0.02 


-3.8 


5.1 


0.13 



Tabic 2. H 2 CS Line Ratios: DR21(OH) 





Ratio 




Value 


T k 
(K) 


N(H 2 ) 

cm -3 


7o,7 

10o,io 

10l,9 


— >■ 6o,6/7 2 ,6/5 ~ > 
— »■ 9o,9/102,8/9 " 
->■ 9l,8/10 3 , 8 /7 " 


62,5/4 
9 2,7/8 
9 3,7/6 


3.93 ± 0.82 
1.15 ± 0.14 
3.42 ± 0.62 


35-33 
23?| 
58l 9 7 


Note a 
Note a 
Note a 


7o,7 — ► 6o,6/10o,io — ► 9o,9 

72,6/5 62,5/4/102,8/9 9 2,7/8 
73,5/4 — ► 6 3i 4/3/10 3i8 /7 — > 9 3)7 / 6 


1.53 ± 0.08 
0.44 ± 0.14 
0.88 ± 0.18 


Note b 
Note b 
Note b 


5 x 10'' 



a. These ratios measure temperature more effectively than density; n(H2)=5 
x 10 7 cm" 3 assumed, b. These ratios measure density more effectively than 
temperature; Tfc=37K assumed.. 

3. Interpretation 

We have used a spherical large velocity gradient (LVG) treatment of radiative 
transfer (RADEX; van der Tak et al. 2007) to model the excitation of H2CS 
within the context of a single temperature and density model of the core region 
of DR21(OH) in which the lines we have observed arise. In this simple model, the 
cloud is constrained to be a very simple spherical entity with uniform spatial 
density, kinetic temperature and molecular column density. The transfer of 
radiation is assumed to be governed by these local conditions. The solutions to 
a statistical equilibrium model of the molecule and the transfer of radiation can 
be expressed in terms of the radiation temperatures of the molecular transitions 
and their associated optical depths. The observed spectrum is then compared 
to the model spectrum. We have included 71 transitions between 34 levels with 
energies below 112 cm" 1 of para-H^CS and 106 transitions between 42 levels of 
ortho-H2CS in our models. 

The collisions of the molecule are a key to modeling statistical equilibrium, 
as collisions provide the link in the model to the local spatial density of collision 
partners, primarily H2. Cross sections for the excitation of thioformaldehyde by 
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collision with molecular hydrogen are unavailable. We have used a scaled version 
of the collisional rates for formaldehyde with He calculated by Green (1991). In 
accord with previous suggestions, we have scaled the rates by a factor of 2.2 to 
approximate the increased interaction with H2 as collision partners. A similar 
treatment has been applied to observations of other massive cores for H2CS by 
van der Tak et al. (2003). 

From the three measured line ratios in the second table, between K=0 
and K=2 lines at 240 and 340 GHz, and between K=l and K=3 lines at 
340 GHz, we find a weighted mean average kinetic temperature of 37±7K. 
From the 7qj — > 6o,6/10o,io $0,9 (E U =91K) excitation ratios, we find a 
mean density of n(H2)=5 x 10 7 cm =3 . The total column density of H2CS is 
N(H2CS)=8.8 x 10 13 cm~ 3 , for the measured linewidth of 5 km s -1 . The ortho 
to para ratio is three. Using these parameters, we calculate the line intensi- 
ties written in Table 1. The agreement with the observations is excellent, with 
the exception of the 10 2) g/8 - ► $2,8/7 hnes, which are too weak by a significant 
factor. The latter are the higher excitation lines in one of the two high excita- 
tion density measures available in our data-the 7 2j 6/5 - ► 6 2i 5/4/10 2j 8/9 — > $2,7/8 
(E U =143K) and 7 3i5/4 -► 63,4/3/103,3/7 9 3 , 7 / 6 (E M =209K) ratios. Unfortu- 
nately, the present cross sections we have employed do not allow a prediction 
of the 10 38 /7 — > 9 3 7/ 6 line. However, we note that while the ratio of the K=0 
lines is high-the lower frequency line is much stronger than the higher frequency 
line-the opposite is true of the higher excitation ratios, where it is the higher 
excitation line which is strongest. The beamsize with which the higher frequency 
lines are measured is smaller than that with which the lower frequency lines was 
measured by a factor of 0.7. If the high excitation emission were to arise in a 
source smaller than the beam, one would expect a relatively stronger intensity 
from the higher frequency transition. 

We conclude that H 2 CS provides a useful tool to measure the temperature of 
warm regions, though some details remain to be explored. CH3CN is often used 
for a temperature probe, but as it is formed in a gas phase chemistry only when 
high temperatures occur, it can be a biased probe. H2CS is formed in an exactly 
analogous fashion to H 2 CO in a chemistry that is fairly temperature-insensitive 
below 300K. It offers advantages over H2CO in the large number of available 
lines, and the much lower abundance (down by more than an order of magni- 
tude). The large number of lines, available to wideband spectrometers such as 
that provided by ALMA, will allow us to deconstruct temperature/density along 
a line of sight. 
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